AN INVESTIGATION OF THE PRODUCTION OF 
MAGNUS FORCES ON A 

HIGH FINENESS RATION BODY OF REVOLUTION 



Timothy Robert Beard 



Library 

Naval Postgraduate School 
Monterey, California 93940 



Monterey, California 





f** 3 * 

Vsasjj. 



AN INVESTIGATION OF THE PRODUCTION OF 
MAGNUS FORCES ON A 

HIGH FINENESS RATIO BODY OF REVOLUTION 

by 

Timothy Robert Beard 

Thesis Advisor: Howard L. Power Jr. 



March 1973 



T 15337 l 



Appnovzd loti pubtic. tidZ^jOit; dliVUbution witimltzd. 



An Investigation of the Production of Magnus Forces 

on a 

High Fineness Ratio Body of Revolution 



by 



Timothy Robert ^eard 
Lieutenant, United States Navy 
B.S., United States Naval Academy, 1966 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 



from the 



% 

> 











Library 

Naval Postgraduate School 
Monterey, California 93940 



ABSTRACT 

A high fineness ratio body of revolution with an ogive 
forebody shape was tested at varying angles of attack, 
Reynolds numbers, and spin rates to measure the Magnus 
forces developed. The experimental data were compared 
with previously derived theoretical equations for an open- 
ended cylinder of revolution to substantiate the theoretical 

work. 
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NOMENCLATURE 



Cy * side force coefficient based on frontal 

D - diameter of missile 

L - length of missile 

p - dimensionless spin rate 

q - freestream dynamic pressure 

R - Reynolds number based on missile length 
Li 

R c - cross flow Reynolds number based on mis 
S - missile maximum cross sectional area 
T - air temperature 
U - freestream velocity 
Y - Magnus side force 
CL - angle of attack 
r - circulation 
fJL. - coefficient of viscosity 
p - air density 
0 J - spin angular velocity 



area 



ile diameter 
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INTRODUCTION 



The study of flight trajectories of spin stabilized 
missiles requires detailed analyses of the aerodynamic 
forces developed during flight. The projectile's lift, 
drag, and side force must be accurately predicted for 
studies of this kind. This paper investigates the flow 
over a high fineness . ratio spinning body of revolution to 
determine the side force and its center of pressure gener- 
ated during flight. The side force causes drift from the 
plane of firing and so is an important aerodynamic force. 

Determining the flow over a missile is a very compli- 
cated three-dimensional problem in which the primary 
mechanism for production of the side force is shed 
vorticity caused by the eventual separation of the boundary 
layer. Because of this the flow problem is dominated by 
viscous effects. Several theories have been developed to 
date to predict the side force and moment on open-ended 
cylinders, assuming laminar boundary layers. These 
theories, however, have yet to be substantiated experi- 
mentally in their appropriate regions of validity. For 
this reason a suitable testing apparatus and test program 
were devised to substantiate these theories. This thesis 
presents work accomplished to date on such studies. 
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HISTORICAL BACKGROUND 



The study of spinning shapes and the side forces 
developed by them has been carried on for many years. To 
properly understand the complicated flow over a high fine- 
ness ratio body of revolution it is useful to study past 
investigations made on less complex spinning bodies. 

In 1671 G. T. Walker [1] recorded the first descrip- 
tion of drift of a spinning body in flight, the curve of 
a sliced tennis ball. Nearly 200 years later G. Magnus [2] 
published the results of his experiments to account for the 
drift of a spinning musket shot. Magnus' method was to 
fire a musket ball whose center of mass was not at its 
geometric center and thus impose a spin on the ball as it 
left the gun barrel. Magnus found a clockwise rotation 
would cause the ball to drift right and a counter-clockwise 
rotation would cause a drift to the left. This substanti- 
ated his theory that an aerodynamic force was responsible 
for the drift. Any such side force produced by a spinning 
body is today called a Magnus force. 

Lord Rayleigh published a paper in 1877 entitled "On 
the Irregular Flight of a Tennis Ball" [3] in which he 
mathematically described the Magnus effects in terms of a 
potential flow around a cylinder with circulation. Missing 
from his theory, however, were skin friction effects 
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between the fluid and cylinder which related the cylinder 
spinning rate to the magnitude of the assumed circulation. 

Little or no theory followed Rayleigh's work until 
after World War I. Shortly after the war a notable attempt 
was made to use the Magnus effect commercially by a German 
named Anton Flettner [4]. He replaced the sails of a ship 
with a large powered rotor. In a cross-wind the Magnus 
effect would power the ship. Two Flettner Rotor Ships were 
built and one sailed across the Atlantic, but the speed and 
reliability of screw propulsion proved to be the rotor 
ship's downfall. 

In 1924 Prandtl [5] published his work on the 
problem of creation of circulation about a spinning body. 

He was able to show that circulation was developed and 
could be accounted for by shed boundary layer vorticity. 

The net circulation was equal but opposite in sign to the 
shed vorticity. 

The Second World War provided graphic evidence of 
the need to accurately predict weapon trajectories. To 
place a shell fired from a gun on target usually required 
more than one shot. Corrections were made from previous 
impacts to achieve the desired results. Significant 
reduction in the number of firings, time, and cost could 
be made if a weapon's trajectory could be accurately 
predicted . 

In 1953 Martin [6] published the first analytical 
method for predicting Magnus forces developed by spinning 
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open-ended cylinders of revolution at small angles of 
attack in incompressible laminar flow. Kelly [7] followed 
in 1954 with a more detailed discussion of Martin's theory. 
Power [8], in 1971, carried Kelly's theory to a high order 
of approximation and included inviscid circulation develop- 
ment to obtain an equation for Magnus force coefficient. 

He improved previous theories by estimating shed vorticity 
that lies close to the body using a cross-flow analogy. 

This theory best models the essential Magnus force contri- 
butions and was used to correlate the wind tunnel data. 

Martin, Kelly, and Power all developed their theories 
by studying open-ended cylinders. Kelly noted in his work 
that there were no known measurements with which he could 
compare his results and that his theory was not directly 
applicable to bodies with different nose shapes. Testing 
open-ended cylinders in a flow is impractical and nearly 
impossible because there is no method to suspend the body 
in the flow and spin it without complicated f low-disturbin 
equipment. For this reason all the testing that is 
described in this paper was done using a cylinder with an 
ogive nose shape. The ogive forebody shape was selected 
because of its constant surface curvature which allows 
theoretical boundary layer calculations to be made more 
easily. 
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MAGNUS FORCE PRODUCTION ON SLENDER BODIES OF REVOLUTION 

AT LOW ANGLES OF ATTACK 



To study three-dimensional flow about a spinning body 
of revolution the full viscous equations of motion must be 
solved. The general equations of motion predict that the 
side force is a function of angle of attack, non-dimensional 
spin rate (the ratio of the surface speed CD r and the free- 
stream velocity U) , body fineness ratio, Reynolds number, 
and Mach number. The complete equations of motion are much 
too difficult to solve, hence simplifications must be made. 
The simplest problem involves incompressible flow with a 
laminar boundary layer, a small angle of attack, and a small 
non-dimensional spin rate. 

A useful method for studying such a flow is by using a 
cross-flow analogy. The basic assumption made to simplify 
the problem is that the flow along the body is independent 
of the flow across the body. This assumption is most 
likely inaccurate for high angles of attack and low body 
fineness ratios. Experimental test data, however, should 
define the region of validity for the application of this 
theory. 

Using the cross-flow analogy, a spinning missile placed 
in a flow at an angle of attack is assumed to have the flow 
velocity about it resolved into two components as shown in 
Figure 1. 
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Figure 1. Velocity Components over a Spinning 
Missile 

The drift experienced by the missile during flight can 
be explained by investigating the effect of the cross-flow 
component (the velocity component perpendicular to the 
missile) . 

From inviscid flow theory it is known that lift per 
unit span on a spinning cylinder in a uniform stream is 

L = p UT . 

where T is the circulation of the flow. The flow pattern 
for the spinning cylinder in inviscid flow is shown in 
Figure 2. 




Figure 2. Streamline Pattern of a Spinning 
Cylinder in Inviscid Flow 
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Inviscid flow theory is helpful for understanding how 
lift is generated on a body in a flow, but it ignores how 
the circulation responsible for the lift is generated. 

Real flow is viscous and does not behave as described 
in inviscid flow theory near the body surface. From bound- 
ary layer theory it is known that the flow around a spinning 
cylinder forms a boundary layer. This boundary layer is 
ultimately shed somewhere on the rear side of the cylinder 
in the flow. Because of the cylinder's rotation, the 
boundary layer formed is not symmetric about the cylinder 
nor is it shed symmetrically in the flow. The boundary layer 
is shed in the form of vortices and is shown schematically 
in Figure 3. 



u r Uq 



Figure 3. Boundary Layer and Shed Vorticity for 
a Spinning Cylinder in Inviscid Flow 

The net vorticity is in the direction opposite to the 
rotation of the cylinder, in this case counter-clockwise. 
Because the total vorticity must be zero for the flow as 
a whole, a circulation in the clockwise direction is gen- 
erated and it is this circulation which provides a 
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major contribution to the Magnus force acting on the 
cylinder . 

In addition to the circulation, boundary layer dis- 
placement thickness, radial pressure gradient in the bound- 
ary layer, and the presence of skin friction shearing 
stresses contribute to the production of a Magnus force. 

The displacement thickness is the distance by which the 
boundary layer would have to be displaced if the entire 
flow were imagined to be frictionless. A representation 
of the displacement thickness is shown in Figure 3. It 
is not symmetric about the cylinder because of the effect 
the cylinder's spin has on the boundary layer growth. 
Knowledge of the displacement thickness permits use of a 
"displaced" body in place of the actual body such that the 
flow around the body can be considered frictionless. As 
can be seen in Figure 3, the magnitudes of the velocity 
components at the boundary layer edge of the advancing sid 
of the cylinder are less than those of the retreating side 
and hence the outer inviscid flow velocity is less below 
the cylinder than above. Because this outer flow is 
frictionless, Bernoulli's equation predicts that the 
pressure below the cylinder is less than that above, and 
a lift is generated. 

Although small, it can be shown that there is a radial 
pressure gradient in the boundary layer around the cylinder 
due to the surface curvature. When integrated, the change 
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in surface pressure due to this pressure gradient contri-* 
butes to the cylinder's lift. 

The presence of friction on the surface of the cylinder 
is also a cause of lift. When the frictional shearing 
stress is integrated over the cylinder's surface there is 
a component of the resultant force in the direction of 
lift. 

Although circulation is a major factor in creating 
lift on the spinning cylinder, the effects of the radial 
pressure gradient and surface friction, however small, 
should be included in the analysis of the total Magnus 
force developed by the cylinder. 



14 



DESCRIPTION OF APPARATUS 



Testing to obtain experimental force data was done in 
the Naval Postgraduate School low speed wind tunnel. The 
tunnel test section was hexagonal in shape with maximum 
dimensions of 3.5 x 5.0 feet. 

The model used for testing was a five- sectioned 
aluminum cylinder 2.950 inches in diameter with an ogive 
nose shape as shown in Figure 4. The missile was con- 
structed in sections so that length- to-diameter ratios of 
9.9, 7.8, and 6.0 could be tested. 

Figure 5 is a line drawing of the model and sting 
balance system. The missile was supported by two Fafnir 
KD66 precision ball bearings which were press-fitted onto 
the sting balance. Each bearing was enclosed in a stainless 
steel collar which allowed the missile to be fastened to 
it \tfith screws. Mounted at the end of the balance was a 
1/8 horsepower Standard Pneumatic rotary-vane, bi-directional 
compressed air motor. Figure 6 shows the balance, bearing, 
and motor configuration. The compressed air was supplied 
to the motor through 1/4-inch copper tubing routed along 
the missile support apparatus. 

Mounted on the drive shaft of the air motor was an 
aluminum drive disk. Located near the disk was an elec- 
tromagnetic tachometer pickup which sensed passage of 
ferrous material (steel screws) in the aluminum disk. 
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Figure 4. Test Model at Positive Angle of Attack 
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Figure 5. Details of Model, Air Motor, and Balance 
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Figure 6. Balance, Bearing and Air Motor 



Also mounted on the balance were two four-element 
strain gage bridge circuits which were used to measure 
strain in the sting balance and thus provide the means of 
determining the Magnus force and its center of pressure. 
Each gage was a Baldwin-Lima-Hamilton strain gage, model 
number FAE 25 12S13, with a gage factor of 2.08. The 
arrangement of the gages for one bridge is shown in 
Figure 7. The gages were wired in such a manner as to 
sense only side force and the Magnus moment applied to the 
balance. Because of the beam geometry and the longitudinal 
arrangement of the gages there was no response to any 
applied torque. Lift forces on the beam also had no effect 
on the bridge output because of gage arrangement. The 
bridge remained balanced when lift forces were applied 
because one pair of gages (#2 and #4 as shown) was in com- 
pression while the other pair (#1 and # 3) was in tension. 
The currents in each pair were different because resist- 
ances were different, but because the voltage drop was the 
same across each pair the bridge remained balanced. 

Because each bridge circuit was made up of four active 
elements the strain indication was magnified by a factor 
of four. The bridge was also temperature compensating 
because all elements were affected the same by changes in 
temperature. Figure 8 is a picture of one side of the 
balance with mounted strain gages. 

The balance was screwed at its base into a tapered, 
circular cross-section, stainless steel sting. The sting 
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Figure 7. Balance Bridge Circuit Configuration 



Tv 
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Figure 8. Balance and Strain Gages 



was held in place at its base by a support and could be 
rotated 90° for measurements in either the horizontal or 
vertical plane. A 90 x 6 x 1-inch aluminum sting support 
was bolted top and bottom to the wind tunnel to support 
the missile and sting arrangement. The beam was machined 
so that the entire sting could be rotated in the vertical 
plane to provide angle of attack as shown in Figure 4. 

The sting support was raised or lowered to keep the model 
centered in the test section. The angle of attack was 
measured with a clinometer. 

Figure 9 is a wiring diagram of the electronics used 

f 

to obtain data. An input voltage of 6 volts was supplied 
to each bridge circuit by a Moxon Electronics SRC Division 
power supply, model 3564. Each bridge output was amplified 
by a factor of 100 by an Aztec model 886 DC amplifier and 
recorded on a Hewlett-Packard two channel chart recorder, 
model 7100B. Two potentiometers per bridge circuit were 
used to initially balance each bridge before testing was 
begun. Figure 10 is a picture of the apparatus used in 
the testing. In addition to the digital counter, a 
General Radio strobe, model 1540-PI, was also used to 
determine missile spin rate to the nearest 10 RPM. 
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Figure 9. Wiring schematic for Electronic Equipment 
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Figure 10. Electronic Equipment 



CALIBRATION AND DATA REDUCTION METHODS 



Airspeed in the tunnel test section was measured using 
water manometry. A pitot tube which measured total and 
static pressures in the test section provided the inputs 
to the water manometry set up. Because the pitot tube was 
located in the most forward part of the test section and 
only six inches off the test section’s floor, calibration 
of the system was necessary. To calibrate the tunnel's 
system another pitot tube was mounted in the center of the 
test section and readings were compared throughout the 
range of airspeeds possible in the tunnel. These corrections 
to tunnel readings were made in the computer data reduction 
program. 

In addition to pitot tube calibration, corrections to 
the test section velocities were necessary due to test sec- 
tion blockage by the model. The correction applied is 
developed in Pankhurst and Holder [9] : 

V true ’ W 1 ' 0 ♦«,> 

where V is the velocity directly measured in the tunnel 
tes 
by: 



test section, is the tunnel blockage factor and is given 



7T i 



= (— ) T (1.0 + . 4D/L) 



Vol 
r 3/2 



7T ^ 

where (-^-) 2 is a factor appropriate to a body of revolu- 
tion, T is a factor dependent on the shape of the tunnel 
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cross section, Vol is the volume of the model, D/L is the 
reciprocal of the model fineness ratio, and C is the tunnel 
test section cross sectional area. 

Corrections to test data due to the effect of air hose 
pressurization on bridge circuit output were necessary. 
Pressure applied to the air motor stiffened the rubber hoses 
leading to the motor and caused a deflection in the sting 
balance which in turn affected strain gage bridge output. 
Corrections for pressurization were made in the computer 
data reduction program to correct for these effects. Input 
data for the pressurization corrections were obtained by 
prior calibration of the sting balance which gave chartline 
deflection for any given air pressure. 

The sting balance system described previously was cali- 
brated by hanging a known weight a known distance from each 
bridge circuit location and recording each bridge's output 
in units of chart lines on the Hewlett-Packard chart 
recorder. Calibration apparatus is shown in Figure 11. 

By plotting the moment versus chartline deflection for 
each bridge circuit linear graphs were obtained whose slopes 
were calculated. Test results (chartline deflections for 
unknown forces) were multiplied by the slopes for each 
circuit to determine moments. By knowing the moment at 
each bridge circuit the applied force and its location as 
shown in Figure 12 were calculated using the following 
simultaneous algebraic equations: 
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Figure 11. Calibration Rig 
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Figure 12. Strain Gage Bridge/Beam Geometry 



F =(M 2 - M 2 )/d 



X = M,/F 
cp 2 
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TESTING PHILOSOPHY AND PROCEDURE 



Testing to determine the Magnus moments was done with 
the objective of obtaining data to verify the theoretical 
equation developed by Power in 1971 [8], which predicts 

the Magnus force coefficient developed by a spinning body 
of revolution. Power's equation showed that the Magnus 
force coefficient was dependent on angle of attack, non- 
dimensional spin rate (ratio of Cylinder surface speed to 
flow velocity), model length- to-diameter ratio, the 
Reynolds number based on missile length, and the cross flow 
Reynolds number. 



R. ' 



: 5 > 2 { 



60.9757 + —h (k [256 C n R 3/4 
r l d 1 C 



357.37] - 227,0112 ( 
R. 



b> ! } 



( 1 ) 



• 6a f 

where C^ = 7Tp ancl given in Figure 12, is a functio 
of cylinder spin rate and angle of attack. By varying 
angle of attack, non-dimensional spin rate, and the model 
length- to-diameter ratio while maintaining constant Reynolds 
numbers for each angle of attack, data could be obtained 
for three length- to-diameter ratios and various spin rates 
which would enable the theoretical equation to be compared 
with experimental data. Because the non-dimensional spin 
rate was in part dependent upon the Reynolds number based 
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Figure 13. Circulation function for spinning cylinders 



on model length, two Reynolds numbers were chosen for testing 
which were compatible with all three model lengths and wind 
tunnel capabilities. The two Reynolds numbers selected were 
700,000 and 1,000,000. 

It was found that pre-determining non-dimensional spin 
rate (p) for the given Reynolds number was impractical 
because the spin rate required for each p was difficult and 
time-consuming to set during testing. Therefore, the pro- 
cedure decided upon was to pre-select various air pressures 
to the air motor and accept whatever spin rates and hence 
p’s that resulted. 

Procedure played a major part in the testing because 
of the time that could be saved through proper planning. 
Because two tunnel velocities and a wide range of non- 
dimensional spin rates (plus and minus) were tested for 
each angle of attack, testing was extremely time consuming. 

By establishing a procedure which involved a minimum number 
of configuration, tunnel velocity, and model spin rate 
changes, the time consumed for any part of the testing was 
reduced to a minimum. 

Upon completion of testing all raw data obtained were 
typed on computer cards and fed into an IBM 360 digital 
computer and reduced by the data reduction program. 
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DATA ANALYSIS 



Due to various physical phenomena encountered during 
the testing the raw data obtained were very difficult to 
analyze to obtain meaningful results. A sample of the 
type of data obtained during the testing from the Hewlett- 
Packard chart recorder is shown in Figure 14. 

The raw data indicated that the missile' vibrated in 
two modes during the testing. The amplitude of the short 
period mode varied between three and five chartlines 
(three and five millivolts) and appeared to be at a con- 
stant frequency throughout the testing. The amplitude of 
the long period mode varied between eight and eleven chart- 
lines in amplitude (measured from the mean of the short 
period oscillation) but its frequency was random. 

Data were obtained by considering the long period mode 
first. The best estimate of its average amplitude was 
made and this information was modified, if necessary, by 
considering the superimposed average short period amplitude. 
In this way a constant voltage value was determined for use 
in the data reduction program. Because of the constant 
frequency of oscillation of the short period mode, averag- 
ing was easy compared to averaging the random oscillations 
of the long period mode. These vibrational modes were 
caused by at least four major effects. First, tunnel wall 
vibration was transmitted to the model through the model 
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support system. The sting and support were bolted directly 
to ' the tunnel floor and ceiling in the test section. Lack 
of time prevented the construction of the originally 
designed isolated support system. Second, mass unbalance 
of the missile itself caused a vibration when the missile 
was rotating. Even though the missile and bearing collars 
were finely machined there was a mass unbalance due to 
varying missile wall thickness and bearing slop. Third, 
asymmetric vortex shedding and leeside separation of the 
flow about the missile added to the vibration. This was a 
function of Reynolds number and rotation rate. Fourth, the 
tunnel test section velocity was not constant and large air 
flow fluctuations occurred. The solution to this problem 
would be to install flow- straightening screens to the tunnel 
to ensure constant low turbulence flow in the test section. 

During the testing a breakdown forced the removal of 
one of the two counter-rotating Bendix constant speed pro- 
pellers used to generate the tunnel's velocity. This 
severely reduced the maximum velocity that could be 
attained and subsequently delayed the testing programs until 
maintenance was completed and a new set of parameters for 
testing, reflecting the change in tunnel speed capability, 
could be computed. The removal of the propeller was nec- 
essary because of lack of parts. The long period of vibra- 
tion was attributed primarily to the tunnel test section 
velocity fluctuations while the short period oscillations 
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were attributed to tunnel wall vibration, mass unbalance 
in the missile, asymmetric vortex shedding, and flow 
separation. 

Because of the magnitude of the oscillations and 

I 

tunnel flow variations, error bounds were computed. Due to 
the type of strip chart paper used (ten chart lines per 
inch) it was assumed that the raw data were accurate to 
within ±1 chartline. It was also estimated that tunnel 
flow fluctuation was ±1 centimeter of water at dynamic 
pressures above 10 centimeters of water and ±.5 centimeter 
below dynamic pressures of 10 centimeters. Evaluation of 
the error bound for each data point was done by computer 
analysis and representative values are shown in Figure 15. 
The error in due to tunnel velocity fluctuations was 
about 111 at each Reynolds number tested. The major error 
contribution was due to the vibration and the resulting 
inability to interpret the strain gage readings accurately. 
The uncertainty bound for each data point i\ras computed by 
adding one chartline to the aft bridge circuit reading and 
subtracting one chartline from the forward bridge circuit 
reading. In this way the maximum testing error could be 
calculated. It can be seen in Figure 15 that the larger 
of these uncertainties occurred for the lowest Reynolds 
number and highest angles of attack (C^ is directly pro- 
portional to angle of attack and inversely proportional 
to tunnel velocity) . Because of the magnitude of the error 
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bounds, interpretation of the data was very difficult. 
Similar problems have plagued all experimenters investigat- 
ing Magnus forces to date. The forces measured are very 
small and the data accuracy required is better than that 
usually required for measuring drag forces. 
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RESULTS 



A representative summary of the data obtained during 
experimentation is presented in tabular form in Appendix B 
and graphically in Figures 16 through 27. Figures 16 
through 20 are plots of corrected data. At each angle of 
attack the data were reduced by subtracting the Magnus 
force coefficient obtained for the no-spin condition from 
each data point. This ensured that the curve of force 
coefficient versus spin rate would pass through the origin 
as predicted by theory. This procedure corrected the data 
for the experimental side force coefficient resulting from 
any missile yawing mis - alignment in the tunnel. Represen- 
tative error bounds for the data points are shown for each 
of the first five plots (16 through 20) . As discussed 
earlier, the bound is larger for the smaller Reynolds 
number and higher angles of attack. Small CL and p theory 
predicts that a plot of versus p is linear. Hence the 
corrected data points in Figures 16 through 20 were fitted 
with a straight line whose slope d^^/dp was cross -plotted 
versus angle of attack for each length- to- diameter ratio 
and Reynolds number. This slope was then adjusted to 
yield a smooth variation for (3Cy/dp versus CL . Justification 
for this method was provided by the size of the error bound 
in each case. 
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Figures 21 through 25 are the resulting plots of the 
Magnus force coefficient Cy versus angle of attack Q . Also 
shown are results developed by Power and Kelly for calcu- 
lating Magnus force coefficients. Figures 26 and 27 are 
plots of Magnus force coefficient versus length- to-diameter 
ratios for constant Reynolds number and angle of attack 
and for various spin rates. 



39 



Re (L) = 7.0x10 



PC Q 

i P, O ^ 

pc :z> 




40 



Figu* 1 . '‘agnus Force Coefficient vs. 
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Figure 17. Magnus Force Coefficient vs. 



Re (L) = 7.0x10 
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Figure 18. Magnus Force Coefficient vs. 



Re CL) = 1.6x10 
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Figure 19. Magnus Force Coefficient vs. 



Re (L) = 1.6x10 
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Figure 20. Magnus Force Coefficient vs. 




Figure 21. Magnus Force Coefficient 
vs. Angle of Attack 
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Figure 22. Magnus Force Coefficient 
vs. Angle of Attack 
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Figure 23. Magnus Force Coefficient 
vs. Angle of Attack 
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Figure 24. Magnus Force Coefficient 
vs. Angle of Attack 
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Figure 25. Magnus Force Coefficient 
vs. Angle of Attack 
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Figure 26. Magnus Force Coefficient vs. L/D 
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Figure 27. Magnus Force Coefficient vs. L/D 



51 



CONCLUSIONS 



The graphical data suggest, but with the error bounds 
of the experimentation cannot verify, that the Magnus 
force might be linear in both 0! and p for small angles of 
attack and small spin rates at constant Reynolds numbers. 
Other experimenters have found non-linear variations of Cy 
at angles of attack greater than 10°. Power's theory pre- 
dicts a linear variation with p and non-linear variation 
with CL while Kelly's [7] predicts a linear variation with 
both p and a . However, due to the magnitude of the error 
associated with each data point obtained, a great deal of 
latitude for data interpretation was available once the 
data were reduced. The data scatter indicates that if the 
sources of error discussed in the data analysis section 
could be reduced, the results would be more valid for 
comparison. 

In general the forces obtained were much lower than 
those predicted by Kelly's or Power's theory as shown by 
Figures 21 through 25. These theories were developed 
assuming an open-ended cylinder while the model had an 
ogive nose shape. This difference alone, however, cannot 
account for the discrepancies between theory and experi- 
mental results. Figures 26 and 27 indicate that as the 
length- to-diameter ratio of the missile increases, so does 
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the side force. However, comparison of the experimental 
data with theory shows that the variation of C y with L/D 
in the region of these experiments is smaller than predicted 
by theory. 

The experimental data presented in Figures 21 through 
25 were used to calculate in equation (1) . The results 
showed to be a small negative number. This implies 
that flow circulation contributed negatively to the Magnus 
force, which contradicts physical reasoning. If C^=0, the 
resulting equation is essentially what Kelly predicted. 
However, the theory must not be discounted simply because 
the experimental data do not support it. When the error 
associated with the data is considered it is readily 
apparent that the data may not correlate with theory. 

Because the measurement of these forces is so critical, 
a better method for measuring them must be devised. 

Although the method used did provide some meaningful data 
which supported theory to some extent, the necessity for 
finer instrumentation and refinement of equipment is 
obvious. Precise machining was necessary for all of the 
apparatus, but no amount of machining or expensive, elab- 
orate electronic equipment will enable accurate data to 
be obtained if other parameters such as poor tunnel flow 
characteristics are a part of the system. The model, as 
constructed, can be tested with a variety of nose shapes 
so that comparisons can be made of the effects different 
forebody shapes have on the production of Magnus forces. 
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However, until more accurate results are obtainable it 
would be useless to test different forebody shapes. Be- 
cause the data obtained do indicate that some side force 
was measured there is sufficient reason for continuing the 
experimentation with re-designed equipment to obtain more 
meaningful results. 
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RECOMMENDATIONS FOR FURTHER STUDY 



1. Design a more sensitive balance that is less 
affected by air motor hose pressurization. 

2. Redesign the sting mounting so that it is iso- 
lated from test section vibrations. 

3. Install flow straightening screens to reduce 
tunnel test section flow fluctuations. 

4. Re-examine the circulation coefficient to 
correlate theory by taking new data with improved exper- 
imental system. 

5. Perform nose shape studies to assess the affects 
of forebody shape on Magnus force. 
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COM PUT FR PROGRAM TO REDUCE RAH DATA TO DETERMINE MAGNUS 
FORCE AND COEFFICIENT WITH ERROR BOUNDS 

DIMENSION ALISTI IOO), BUST! 100). FUST (100) tGL 1ST < 100) » 
PLIST(IOO) , CL I ST { 1 00 ) 



EXPLANATION OF 
XL/XLF - 
DF 
XLD 
T/TR 
P 

RHO 

XMU 

R 

G 

XMF/XMA - 
DF 

AOA/ARAD - 
0 

V/VTRU - 
EPS 
TSEC 
XCO 
CONE 
VOL 
REV/XO 

REL/REC - 

XCO 

S 

FCOR/ACOR - 



SYMBOLS USFD IN T H E PROGRAM 
MISSILE LENGTH ( IN. /FT. ) 

MISSILE DIAMETER IN FEET 
MISSILE L ENGTH-TO— DIAMETER RATIO 
TEMPERATURE ( FA HP ENHE IT/RANK I NE ) 
ATMOSPHERIC PRESSURE (IN. HG.) 
ATMOSPHERIC DENSITY 
COEFFICIENT OF VISCOSITY (AIR) 

GAS CONSTANT 
GRAVITY CONSTANT 

SLOPES CF MOMENT VS CHART DEFLECTION 
FROM STING CALIBRATION 
DIAMETER CF MISSILE (FT.) 

ANGLE OF ATTACK ( DEGREE 3/ R AD I ANS ) 
DYNAMIC PRESSURE IN CM H2G 
TUNNEL VELOCITY ( ME ASUP.ED/CORRECTED ) 
VELGCITY CORRECTION FACTOR 
TEST SECTION CROSS SECTIONAL AREA 
LENGTH OF MISSILE NOSE SHAPE 
VOLUME OF MISSILE NOSE PIECE 
VOLUME OF MISSILE 

REVOLUTIONS PER MINUTE/RADIANS PER 
SECOND 

REYNOLDS NUMBERS BASED ON LENGTH/CROSS 
FLOW 

LENGTH OF MISSILE NOSE SHAPE 
MISSILE CROSS SECTIONAL AREA 
CORRECTIONS TO DATA DUE TO EFFECTS OF 
HOSE STIFFENING ON STING 



WRITE (6, 10) 



R F AD IN CIRCULATION FUNCTION DATA FOR USE IN POWER'S 
THEORETICAL EQUATION FOR MAGNUS FORCE COEFFICIENT 

READ( 5. 14) (PL1 ST(J ) , J = I ♦ 1 7 ) 

RE AD (5. 14) (CLIST(J ) » J = 1 , 1 7 ) 

DO 100 1=1,10 
READ(5.U)L,NDATE.M,E 
IF (E.EO. 1.0JG3 TO 9 
IF(E.GT. 5. 0)GD TO 1 
WRIT? (6, 12 ) NDA T E 
GO TO 2 

1 WR IT F ( 6 , 13 1 NDA T E 

RF AD IN DATA FOR FORCE READING CORRECTIONS DUE TO HOSE 
AIR PRESSURES 

2 RFAO( 5 . 14) ( ALI ST( J) , J = 1 , L) 

RFAD15, 14) (BLIST (J ) ,J=1,L) 

RPAD(5 . 14 ) ( FLI ST ( J ) ,J = 1,L) 

R F A D ( 5.14) (GLI ST £ J ) ,J=1,L) 



DO 101 N =1 . M 

READ IN MISSILE LENGTH. TEMPERATURE (DEGREES 
FAHRFNHEIT), ATMOSPHERIC PRESSURE, AND A MISSILE 
PARAMETER USED TO DETERMINE LOCATION OF CENTER OF 
PRESSURE 

READ ( 5. 15) XL.T.P.D.K 
XL F = XL/12. 
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DF = 2.950/12. 

XLD = XL F /DF 
WRITE(6. 16 )XL, XLD 

TR = T +459. 6 7 
R = 53.34 
G — 3 2 \.~~l / 1 

RHO = Pv-2116. 22/( 29. 92"' G* R*TR) 
X = 2.26967 F-08 
XMU = X? ( T R L L 1 . 5 ) / ( T R + 1 9 8 . 7 2 ) 
XMF = 63.7/274.8 
XMA = 86.4/323.7 



CD 102 J = 1 , K 

R F AD IN MISSILE ANGLE OF ATTACK, DYNAMIC PRESSURE IN 
CM H20 (CORRECTED TO TRUE 0 IN PROGRAM), THE 
REVOLUTIONS PER MINUTE, AND THE TWO CHARTLINE 
READINGS (1 CHARTLINE PER MILLIVOLT) 

READ( 5, 17) AOA, 0, REV , VF , VA, PHOS E 
SPIN = REV 
ANGL = A8S ( AOA ) 

ARAD = ANGL'’ 3. 141 59/180. 

CA = COS(ARAD) 

SA = SIN(ARAD) 

Z = 0*62.441 . 13 3 / ( RHO? 15 . 24 ) 

V = SORT ( Z) 



INSFRT DATA TO DETERMINE VELOCITY ERROR BOUNC 

IFIO.GT.IO .00) GO TO 3 
OX = 0+0.5 
GO TO 4 

3 OX = 0 + 1 .0 

4 ZP = 0X462.4-1 . 123/( RHOL 15.24) 

VP = SORT ( ZP ) 

TSEC = 14.41 
CONE = .01973 

XCO = 7.244/12. . '■ 

VOL = CONF+3. 14159*(DF**2)* (XLF-XCOJ/4. 

EPS = .75--( (3. 14159/4. )*'>: . 5) s ' ( l. + .4vDF/XLF)4V0L/ 

(TSEC?* 1.5) 

VTRU = V4( l. + EPS) 

VTRP = VP? ( l.+ EPS) 

5 XO = REV?- 2.4 3.14159/ 60 . 

6 PB = DF4X0/VTRU 

P BP = DPtXO/VTRP 

RFL =• RHG? XLF4 V7RU+CA/XMU 
REC = RHG* DF? VTRU? SA/XMU 
' S = ( (DF/2. )' 42 ):-3. 14159 

FCOR = PIF2( PHOSE, A LI S T , L , B L I ST ) 

VFC = VF-FCGR 

ACOR = PIF2(PH0SE ,FLI ST ,L,GLIST) 

VAC = VA-ACCR 
XMOF = XMF? VFC 
XMGA = XMA? VAC 



DATA INPUT FOR DETERMINING ERROR IN INITIAL RAW DATA 
READINGS 

VF P = VFC + 1.0 
VFM = VFC- 1.0 
VAP = VAC + 1.0 
VAM = VAC-1.0 
XMFP = XMF4VFP 
XMFM = XMF? VFM 
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XMAP = XMA* VAP ' 
XMAM = XMA*VAM 



CALCULATION OF MAGNUS FOPXE, FORCE COEFFICIENT, AND 
ERROR BOUND 

Y = ( XMOA-XMCF) /5 . 

XCP = XMOF/Y 

CY = 2 .*Y/(RHOT VTRU-WTRUOS) 

YP = ( XMAM-XMFP) /5. 

CYP = 2.VYP/C RHC^' VTRU" VTRU :i? S) 

YM = ( XMAP-XMFMJ/5. 

CYM = 2.^YM/( RHC* VTRU“ VTRUTS) 

CYVP = 2.--Y/(RH0 ! -VTRP^ VTRP« S) 

IF(AOA.LT.1.0)GO TO 7 
if(pb.f:q.o.o)go to i 

PSIN = A 3S ( PB/ S A } 

F = PIF2(PSIN,PLIST,17,CLIST) 

Cl = .6-‘ APAD--F/( 128. *3. 14159? PB) 

CYP W = ( ARADC-PBTXLCT XL D/ ( R EL** .5 ) ) S ( 60 . 97 + ( 256 ,*C 1$ 

( RFC*? . 75 3-357.37 )*XLD/( REL* t, 5 )-227. 01^ XLDAXLD/REL) 
GO TO 8 
7 CYP W = 0.0 



PRINTOUT DATA DESIRED 

8 WRITF (6 , 18 ) AOA ,REL, VTRU, SP I N , P B , Y, C Y ,C YP , C YM , C YV P,C YPW 

102 CONTINUE 
101 CONTINUE 
100 CONTI Nil F 



FORMAT CONTROL STATEMENTS ARE: 

10 FORMAT! ///T40, ’APPENDIX B' ///T31 MAGNUS FORCE AND', 

IX, ' COEFFICIENT • ) 

11 FORMAT ( 31 5. F5. 1 ) 

12 FORMAT (////AX, 'TESTING DA TE : ' , I 3 , 1 X , • JA NUA RY 1973') 

13 FORMAT (////AX, 'TESTING DAT E : • , I 3 , 1 X , • FE B RU ARY 1973') 

1A FORMAT (7 FI 0.1) 

15 FORMAT (F10. 3, F10. 1,2F 10.2,15) 

16 F0RMAT(//5X, 'MISSILE LENGTH = ' , F6 . 3, IX, ' I NCHES • , 1 X, 

•(1/0 = ' ,F3. 1 , ' ) ' //6X, ' AOA* , 6X, ' RE(L) ' ,7X, • V' ,5X, 

’ RPM* , 3X, • P-HAT ',3X,'Y',5X,'CY',3X,'CY(-)',2X,'CY(+)' 

. 2 X , ' CY ( V+ ) • , IX, ' CYPOW ) 

17 FORMAT ( F10.1,F10.2,AF10. 1) 

18 FORMAT (' O' , AX, FA. 1, IX, IP E 1 2 . 5 , 2X , OP F 5 . 1 , 2X , F6 . 0 , 2X , 

FA. 1,2X,FA.2,7(2X,F5.2) ) 

19 FORMAT ( 7F10.3) 

9 STOP 
END 



FUNCTION PIF2(X, FL 1ST, N, GUST) 
SUBPROGRAM FUNCTION 2ND ORDER LOCKUP 
DIME NS I CN FLIST( 100) *GLIST (100) 
BLIF(P,Q,R,S,T ) = (( 0-P ) ( S-T ) / ( R-0 ) +S ) 

IF (X-FLIST(N) ) 2,1,1 

1 I = N— 1 
GO TO 5 

2 I F( X-FLI ST ( 1 ) ) A, A, 6 
A 1=1 

5 K = 1 

GO TO 30 

6 K= 2 
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7 

8 

9 

30 

10 

11 

12 

13 

14 

15 

16 

17 

18 



DP. 8 I = 1 , N 
IF (X-FLISTU) ) 9, 9,8 
CD NT INUE 
I =N 
1 = 1-1 

BLIF1=BLIF(X,FLIST(I) .FLISKI + 1) ♦GLIST(I ) .GLISKI + l) > 
IF ( K-l ) 11,11,12 

PI F 2 = BLIF1 
RETURN 

I F ( d +2 )— N ) 13,13,16 

I F ( (I — 1) — 1) 15,14,14 

IF (ABS ( FLIST( I-l)-X) -ABS ( F L I ST (I + 2 ) -X ) ) 16,15,15 

L = 1+2 
GO TO 17 
L= 1-1 

BLIF2=ELIF { X.FLIST ( I ) , FlIST ( L ) , GL 1ST (1) , GL 1ST (L ) ) 

P I F2=BL IF(X,FLIST(I + 1) ,FLIST(L) ,BL1F1,BLIF2J 

RETURN 

END 



I 
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